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Abstract. We derive an exact analytic solution to a Klein-Gordon equation for 
a step potential barrier with cutoff plane wave initial conditions, in order to 
explore wave evolution in a classical forbidden region. We find that the relativistic 
solution rapidly evanesces within a depth 2xp inside the potential, where Xp is the 
penetration length of the stationary solution. Beyond the characteristic distance 
2xp, a Sommerfeld-type precursor travels along the potential at the speed of light, c. 
However, no spatial propagation of a main wavefront along the structure is observed. 
We also find a non-causal time evolution of the wavefront peak. The effect is only 
an apparent violation of Einstein causality. 
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1. Introduction 

Since the early beginnings of quantum mechanics, the problem of particle propagation 
in classical forbidden regions has been the subject of both theoretical and experimental 
investigations. Over the years, several non-relativistic approaches based in cutoff 
wave initial conditions have been introduced in the literature in order to investigate 
the time-dependent features of wave evolution in evanescent media. Some of these 
theoretical models H, ^ were inspired in the pioneering work of Sommerfeld 
and BriUouin [ |6| % while others |9[ |lO[ ^ were based in the seminal work of 
Moshinsky p3|, |14[, who a few decades ago started a fundamental discussion on the 
non-relativistic and relativistic transient effects. These models represent important 
steps towards the clarification of the dynamics in classical forbidden regions, and 
a renewed motivation to explore this problem has been recently stimulated by the 
issue of superluminal velocities in photon ^ and microwave 

HI tunnehng. 

Hence, it is clear that a full relativistic approach to describe the wave evolution in 
evanescent media is needed. Nevertheless, this has become a complex problem due to 
the lack of exact analytical solutions to relativistic wave equations with appropriate 
initial conditions. Among the few works in the field ^ Deutch and Low 
have provided a lucid description of barrier penetration of an initial state given by 
a cutoff Gaussian wavepacket, based on a one-dimensional Klein-Gordon equation. 
Although no exact relativistic solutions were obtained, the issues of Einstein causality 
and superluminal phenomena were rigorously discussed using approximate solutions. 

In this paper we consider a model based on the Klein-Gordon equation, as in the 
work of Deutch and Low , for a potential step barrier and cutoff plane wave initial 
conditions. We obtain an exact analytic solution to the problem along the potential 
region and study the main features of wave evolution, in particular the regime of 
transient effects at early times. 
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The paper is organized as follows. In section 2 we discuss the shutter problem, 
and present the analytical derivation for the solution to a relativistic wave equation 
for a step potential barrier. Section 3 deals with a numerical example for the solution 
along the internal region of the potential, and the results are discussed in section 4. 
Finally, in section 5 we present the summary and conclusions. 



2. The relativistic shutter problem 

To investigate the time evolution of cutoff plane wave in a classical forbidden region, 
let us consider a classical field ip^ satisfying a one-dimensional Klein-Gordon equation 
with a variable potential V{x), as in the model of Deutch and Low 
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-^1pl{x,kr,t) = -^-^lPl{x,kr,t) + Vo{x)if^{x,kr,t). (1) 

In our case, Vq{x) is given by a step potential barrier, 

III 

where /zq = (moc/h), and the initial condition at i = corresponds to a plane wave 
shutter (see figure ^, given by. 

The simplicity of our quasi-monochromatic initial state (||), allows a closed 
analytical solution of the problem. It differs from that of reference where a 
cutoff Gaussian wavepacket initial condition was considered. Note that condition (||) 
comes from the fact that for i < 0, the solution for the left side of the shutter is 
given by iprix, kr,t) = eKp[ikr{x — ct)] — CKp[—ikr{x + ct)], for a; < 0, and zero for 
a; > 0. 

To obtain the solution for a; > and t > 0, we shall proceed along the same lines 
as in our recent work pT[ |. We begin by Laplace transforming the equation using 
the standard definition 

/>OC 

^{x,kr,s)^ / i}}{x,kr,t)e~''^dt, (4) 

with the initial condition given by equation (^. As a consequence, one gets a pair 
of differential equations corresponding to the regions a; > and a; < 0. In order to 
obtain the transmitted wave function, one must consider the matching conditions for 
the wave function and its derivative at x = 0. The Laplace transformed solution reads, 

?:(^, s) = ^ , e-^-/-, (5) 

I (s + lE) [s+p) 



where p = {s^ + p^c^Y/"^ , and E — kr — (Er/hc) corresponds to the relativistic energy 
Er given in reciprocal units of length. 
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The time dependent solution for a; > is readily obtained by performing the 
inverse Laplace transform of equation (|^) using the Bromwich integral formula, 

7'+'too 

V';(^,i) = 2^ / (6) 

where the integration path is taken along a straight line s = j' parallel to the imaginary 

axis in the complex s— plane. The real parameter 7' can be chosen arbitrarily as long 

as all singularities remain to the left-hand side of s = 7'. 

The integral (|^) expressed in this form, is difficult to manipulate since the integrand 
has branch points at s = ±j/ioc. To surmount this difficulty, let us introduce the 

change of variable, — m = (s +p)/fioc, which allows to eliminate the branch points. 

Thus, p = i^Qc{u^^ — u)/2, and as a consequence, the Bromwich integral may be 

written as, 

27+00 

i^ix,t) = ^ f F{u)du, (7) 

2711 J 

i'y — oo 

where the new integrand F{u) is given by, 



F{u) 



2E jl-u^) 

fiQ u^{u'^ — 2Eu/fio + 1) 

X exp{iiio[u{x — ct) — u^^{x + ct)]/2}. 



(8) 



Note that the branch points go into an essential singularity at u = and two 
simple poles u± — {E i: iq)/ fj,Q, where we defined q = {^0^ ~ E"^)^/"^ . The integration 
in equation (|^) is performed along a straight line L parallel to the real axis cutting the 
positive imaginary axis at 17. We proceed to evaluate the above integral by considering 
a closed Bromwich integration contour (see figure ||), and Cauchy's residue theorem. 
For a:; > ci we close the integration path L from above, by a large semicircle Fi of 
radius R, forming a closed contour Ci. The contribution along r2 vanishes as i? ^ 00, 
and since there are no poles enclosed inside Ci, tp{x,t) = for x > ct. For the case 
X < ct, we close the integration path from below with a large semicircle F2. The 
closed contour C2 contains three small circles Co, C+ and C_ enclosing the essential 
singularity at w = and the simple poles at m+ and u_ , respectively. Hence by using 
Cauchy's theorem, it follows that. 
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2TTi 



17 + 00 



Co C-f. 



F{u)du = 0. 



(9) 



The integrals corresponding to the contours C+ C_ can be easily evaluated, and 
yield the exponential contributions to (|^), namely. 



/ F(u)du = k±e 

2m J ^ ' 



{^qx — iEct) 



(10) 



c± 
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where we defined k± — 2E/{E ± iq). 

The contour integration for Cq requires a more elaborate calculation, since involves 
an essential singularity at it = 0. For this case, we introduce the change of variable 
given by w = — m^^^, thus the integral now becomes 

F{u)du= / — — ^ duj, (11) 

Co C', 

where uj± — [E ± iq)/ifj,o^- To carry out the integration, first let us separate the 
integrand into partial fractions, and substitute the well known formula for the Bessel 
generating function. 



g,(<.-^-)/2 ^ ^^n^^(^) ^ ^(-1)"C.-"J„(,7) (12) 
n— n— 1 

and the series expansion, 

{lo± - Lo)-' = (c.±)-i Y.iuj/LU±rj^{r^). (13) 

n=0 

The resulting integrals can be evaluated by means of the residue theorem. For the 
case of an essential singularity, the residue may be determined by computing explicitly 
the coefficient corresponding to lo^^ from the series expansion and their products. In 
that case, equation ([III) becomes. 



1 
2^ 



F{u)du : 



Co 



2iE 



Ji(r;)-fc+y(-l)"^^ 



n=ll ^ ' 



(14) 



Finally, substituting the results given by equations (14) and ( p^ into equation 
, the solution for the internal region is. 



with tp±{q) defined as, 



i^+{q) + ij-{q), t>x/c 

0, t < x/c, 



(15) 



n=0 



(16) 



In the above expression, Jniv) stands for the Bessel function of order n. The other 
parameters are defined as: ^ = [{ct + x)/{ct — x)]^^'^ , r] = f^o^c^t"^ — x^)^!'^ and 



z± — (E' ± iq)/^Q. From equation ([15[) we see that the solution obeys Einstein 
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causality, i.e. no propagation faster than the speed of hght, c, is detected along the 
barrier region. In other words, an observer located at an arbitrary position xq inside 
the barrier must wait a time t = (xo/c) before detecting the arrival of the signal. 

For the sake of completeness, let us now consider the asymptotic behavior of 
■4>'l{x,t) for the cases /xq ^ 0, < — > oo and x ct. From the solution we have just 
discussed, one may recover the free propagation solution in the limit /zo — > 0. This 
corresponds to let the variables rj 0, q lE. To illustrate the limit process in 
equation ( |l5|), let us rewrite the solution by using equation (p^), namely, 



rrM = k^ [e('^-'^^*) - Jo(77) 



OC 

E 



(CAz-)"J„(??) 



^(z+/ze)"J„(77) 



(17) 



As ^0 — * 0, the variable Jo(?7) — > 1, and since (z_)~^ the first series on the 
right hand-side clearly vanishes. It can be shown that the second series also vanishes, 
by replacing the Bessel functions by their asymptotic values for small values of the 
argument 77, 



J„(,7) ~ 2-"r;"/n!. 
Therefore, one obtains the solution for the free propagation case, 

^ikr{x — ct) 1 



i)l{x,t) 



0, 



t > x/c, 
t < x/c. 



(18) 
(19) 

(x/c) fulfilling 



Note that the free case solution rises from zero only after a time t ■ 
relativistic causality, and then oscillates periodically thereafter. 

In the case of the long-time limit (t — > 00), we have ^ — > 1 and rj — s- 00. From the 
asymptotic expansion of Jn(jl) for large values of the argument rj, 



(20) 



and therefore Jn{v) ~* 0- O^i^ can see from equation (16) that the series in ip+{q) 
vanishes. On the other hand, if we rewrite ?A-(g) by means of equation (p^), 
the exponential term is canceled and the remaining series vanishes. Consequently, 
ij}^{x,t) goes into the stationary solution (/)f.{x,t) given by. 



(21) 



The asymptotic behavior near the relativistic cutoff, is obtained when x ^ ct 
in equation (|l|). In this case we have 77 ^ 0, which allows us to substitute the 
asymptotic expansion ( p8| ) in '>p+{q) (equation (16)). Thus, the series of equation (16) 
goes into an exponential function, which cancels exactly with the exponential term, 
and as a result the solution ip+il) goes like iEJi{ri)/fio^. From similar considerations 
on •il}-{q), an identical expression is obtained and the approximate behavior of ipr{x, t) 
near the cutoff is given by. 



Mo? 



(22) 
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where for exactly the value x = ct, the above expression goes to zero fulfilling 
relativistic causality. 

3. Examples 

In order to exemplify the evolution of the solution given by equation ( p^ ) along the 
evanescent region, we decided to study the properties of \tp^{x,t)\^ as a function of 
time t and the position x. The parameters for the system considered in all the cases 
for the present study are: barrier height /xq = 1.542 nm~^ , incidence energy Ej. — 10.0 
eV {E^ 5.064 x lO'^ nm'^). 

The first case corresponds to the spatial evolution of IV'r {x, t)\^ along the dispersive 
region. In figure |^ we show at early times the birth of the main wavefront as a function 
of the position, for increasing values of time: ti — .001 /s, t2 — .0035 fs and ^3 = .0075 
fs. The solution rises as time goes on, and at ^3 , \^^{x,t)\^ has already crossed over the 
stationary solution [ (a; , f ) ^ (dashed line) . The inset of figure || shows the crossover of 
at later time ^4 — .012 fs. This behavior is relevant since it indicates that 
the relativistic solution fluctuates around |(/>^(a;,t)p, before reaching its asymptotic 
regime. At the inset, we can also observe how the solution evanesces within a finite 
depth given approximately by 2xp = 1.317 nm, where Xp = (I/9) is the penetration 
length of the stationary solution — |fc+|e~*^ (equation (pT])). We find that 

beyond 2xp the solution exhibits a small maxima, corresponding to the birth of a 
forerunner. In figure ^ we depict the spatial evolution of |^^(a;,t)|^ (solid line) for a 
fixed value oit — 0.05 J s. As can be seen in this example, the main part of the wave 
rapidly evanesces in the potential region for small values of the position. However, 
from approximately 2xp onwards, the solution exhibits an oscillatory behavior before 
reaching the relativistic cutoff at a; = 15.0 nm, corresponding to the earliest arrival of 
the signal at a point located within the potential. The stationary solution |</)*(a;,t)p 
(dashed line) is also included for comparison. It is interesting to note the similarity of 
the oscillatory structure in figure ^ to the well known Sommerfeld precursor |2l[ |, 
which is one of the essential features of wave propagation in dispersive media. Despite 
the fact that Sommerfeld's approach is quite different from ours, the similarities go 
beyond the numerical results. For instance, their asymptotic analysis showed that 
the wave function is governed by a first order Bessel function near the relativistic 
cutoff. Our analysis reproduces such behavior, which is given by equation (^). 
For comparison, the value of the Bessel function Ji(?7) modulated by the prefactor 
2iE/ fiQ^, is also included in figure^ (dotted line). Note that if we define the frequency 
of the oscillations of the precursor in terms of the distance between successive zeros of 
Ji{r]), one sees from the definition of the argument i] that the value of the frequency 
depends only on the position x and the potential ^0 that characterizes the medium 
i.e. the precursor frequency is independent of the incidence energy. 

In figure ^ we show \ipf.{x,t)\'^ (solid fine) as a function of the position x, at a 
subsequent time t = 0.3 fs. We can see that the solution reaches its stationary value 
\(j)^{x,t)\'^ (dashed line) for small values of x; nevertheless, near the relativistic cutoff 
at X = 90.0 nm, the precursor exhibits a rich oscillatory structure. The inset of figure 
^ illustrates the forerunner near the cutoff, and shows that the asymptotic behavior 



is dictated by the Bessel function of equation (22) (dotted fine). 

Up to here we have illustrated the spatial behavior of \'ip!^{x,t)\'^ , and some 
interesting features of the time evolution. In order to fully explore the relevant features 
of the time evolution, in figure || we plot \ipri^^ a function of time at different 
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positions: xi — 0.4 nm , X2 — 0.6 nm and X3 = 0.8 nm. For all the curves depicted, 
as soon as i > (x/c) the solution is different from zero along the internal region, 
fulfilling relativistic causality. As can be seen, the solution rises from zero at t = {x/c) 
and grows monotonically towards a maximum value, from which it starts to oscillate 
thereafter, forming a pattern very similar to the diffraction in time phenomenon p^ . 
The concept of diffraction in time was originally introduced by Moshinsky [|l3j while 
discussing the shutter problem for the free particle Schrodinger equation. He observed 
a time-dependent oscillatory regime of the probability density near the semiclassical 
wavefront that he named diffraction in time, in analogy to the well known Fresnel 
optical diffraction. It is interesting to note the resemblance of the oscillatory pattern 
in figure ^, to the diffraction in time phenomenon observed in the free propagation 
11 . Moreover, in the low-energy regime {^.q/E ^ 1) the solution ( |l5| ) can be 



rewritten in a more concise form by using equation ([l2|), namely, 

« 2{E/V) [C/3(*??/^,7?) - U,{^v/^,v)] , (23) 

where Ui and C/3 are the Lommel functions of two variables p^ , widely used in 
connection with optical diffraction p^ . The resemblance to diffraction phenomena 
suggests that there exists a more profound link; however, the physical implications of 
the striking mathematical similarities found above deserves further study. 

It is important to mention that the transient effect depicted in figure |6| is observed 
in the low-energy regime i.e. {no/E) ^ 1; this condition is satisfied in the present 
example, where the effect was observed for values of the ratio {no/E) ~ 30. Moreover, 
we only observed the phenomenon in the regime of small values of the position x where 
the solution decays in the potential region i.e. x < 2xp. From values greater than 
X ~ 2xp the solution enters into a different oscillatory regime, and the diffraction-type 
pattern begins to disappear. In figure |^ we illustrate the inhibition of the diffraction- 
type pattern for a fixed value of the position x = 3.0 nm. Clearly, it fades out and is 
replaced by a series of oscillations, which register the fast crossing of the precursor at 
X = 3.0 nm, and the remaining forerunners. 

There is another interesting feature in the time evolution of \ip!^{x, t)p that can be 
appreciated in figure^, in which we plot \'ijj^{x,t)\'^ as a function of time in the main 
peak region. Surprisingly, the maximum peak of the wave appears on x = 0.5 nm 
(dotted line) earlier than the peak at x = 0.3 nm (dashed line) and x = 0.1 nm (solid 
line) . This relative time shift of the wave peak is an apparent violation of relativistic 
causality, and can be interpreted as a non-causal behavior. This comes from the fact 
that we are comparing the maximum wave peak at different positions. However, we 
observe that the wavefront always fulfills Einstein causality, and no signal travels faster 
than c in the dispersive region. Therefore, the observed shift of the main peak may 
be interpreted as a reshaping of the wave and not as a genuine violation of relativity. 

It is interesting to mention that we have observed a similar non-causal behavior in 
the probability density along a classical forbidden region of a rectangular potential 
barrier, within a non-relativistic framework. Moreover, some authors have also 
reported non-causal phenomena in electromagnetic evanescent modes |]l8|. 



4. Discussion 



The possibility of describing the wave evolution from the transient to the stationary 
regime, offers a clear advantage over the asymptotic methods of solution available in 
the literature, for which the short and intermediate transient regimes are inaccessible. 
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In what follows, we shall discuss the new features in the dynamical process of 
evanescent waves observed in the previous section. 

The buildup of ji/'rP exhibits a very interesting behavior; the solution instead of 
just grow monotonically towards the stationary solution I^^Pi fluctuates around such 
value before reaching the asymptotic regime, as it is shown by the series of curves 
of I'f/'^P versus x at different times (see figure ||). The effect of these fluctuations in 
a plot of IV'rP versus t {x fixed) is manifested as a series of oscillations similar to a 
diffraction in time pattern, see figure ^. The inset of figure || shows that beyond a 
certain distance, an interesting structure of the wave appears; this is the birth of the 
Sommerfeld-type forerunner which travels along the potential region, as illustrated 
in figures ^ and ||; the head of this signal propagates at the speed of light and can 
reproduced by the first order Bessel function Jx{r]) (dotted line). The birth of the 
forerunner is an important event since its propagation at longer times becomes the 
dominant process; this is also the case in the context of different relativistic and non- 
relativistic approaches |^], where the characterization of the forerunners has been 
recently emphasized. 

At early times and for small values of the position, the main front of the 
wave decays exponentially along the potential. As time goes on, the dynamics is 
dominated by the propagation of the forerunners since the main front rapidly reaches 
its asymptotic value without propagation. Thus, one may speak of two regimes, 
which as discussed in the previous section, are characterized by 2xp where Xp is the 
penetration length. If wc choose a position x > 2xp and wait for the main wavefront, 
instead of detecting its arrival we would only register the fast crossing of the precursor 
(see figure ^ . The absence of main wavefront propagation in the evanescent region is in 
agreement with a series of works |2^, ^ , which have questioned the existence of 
semiclassical wavefront propagation proposed by Stevens and supported by Moretti 

i'i- 

Another important result of this work is the non-causal peakshift exhibited in 
figure ^ The non-causal behavior observed here is a consequence of the reshaping 
of the wave; reshaping effects have also been observed in the context of wavepacket 
evolution within both relativistic ||l9| and non- relativistic ||2^ approaches. The role 
of this effects and the issue of non-causal behavior has been recently discussed by 
Deutch and Low [|9| for the case of Gaussian wavepacket evolution in the transmitted 
region of a potential barrier, based on approximate solutions to the Klein-Gordon 
equation. Although the barrier and the step potential are different systems, both 
exhibit an evanescent region; hence, we believe that the non-causal behavior observed 
in reference could be related to a reshaping process occurring inside the barrier, 
similar to the reshaping observed in the step discussed in our model. However, in order 
to investigate such a reshaping inside the barrier, the solution of the Klein-Gordon 
equation is required for the internal region. In this respect, the analytical techniques 
used in this work may provide a suitable method of solution to tackle this fundamental 
problem; nevertheless, this not an easy task since the extension of our model to the 
case of a barrier of finite width involves more complicated analytical properties of the 
solution due to the presence of resonances. 

5. Summary and conclusions 

We have derived an exact analytical solution to a Klein-Gordon equation for a step 
potential barrier, using a cutoff plane wave initial condition. To our knowledge, this is 
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the first model which allows a closed solution for the description of relativistic transient 
effects in a classical forbidden region. 

The main features of the spatial and time evolution along the evanescent region 
can be summarized in the following points: (i) We found a regime where the solution 
is exponentially suppressed and thus, decays as a function of x along the potential. 
This main part of the wave does not propagate along the structure. The regime is 
characterized by a region extending from a; = to cc ~ 2a;p, where Xp — (1/g) is the 
penetration length of the stationary solution However, from onwards, 

the solution exhibits an oscillating pattern near the relativistic cutoff, traveling at the 
speed of light, c which can be identified as a Sommerfeld-type precursor. Also, within 
the finite depth 2xp, we found in the low-energy situation that the time evolution of 
exhibits a transient effect similar to the diffraction in time phenomenon 
|l3| . (ii) We showed that along the internal region, there exists a time shift associated 
to the main peak of the wave function, that can be interpreted as a non-causal behavior 
along the classical forbidden region. This of course is only an apparent violation of 
relativistic causality, since in our model the wavefront satisfy always Einstein causality 
i. e. no signal travels faster than the speed of light. 

The relevance of these results, comes from the fact that our findings may be of 
interest to elucidate on the problem of wave propagation in finite width potentials. 
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Figure captions 



Figure 1. Shutter problem for a potential step barrier Vq. An initial state ^l>{x,0} 
in the region a; < is instantaneously released at t = by the removal of the shutter 
S. 



Figure 2. Integration contours C2 = L + r2 + Co + C+ + C— and Ci = L + Fi , used 
to evaluate Eq. (7). The infinite semicircles Fi and F2 (dashed line) correspond to 
the cases x > ct and x < ct, respectively. 



Figure 3. The birth of \ip"(x, t) p (solid line) as a function of distance x for increasing 
values of time: ti =0.001 fs,t2 = 0.0035 /s and = 0.0075 /s. Note that ji/"^ (x, t)|2 
fluctuates around the stationary solution \<f>f,(x,t)\'^ (dashed line). The inset shows 
at a later time ti = 0.012 fs the birth of a Sommerfeld-type precursor near the 
relativistic cutoff at a; = 3.0 nm. 



Figure 4. Plot of \iIj"{x, t)p (solid line) as a function of distance x for a fixed value 
of time t = 0.005 fs. Notice that the wave function exhibits a Sommerfeld-type 
precursor near the relativistic cutoff at x = 15.0 nm. The precursor is accurately 
described in the vicinity of x=ct by a Bessel function (dotted line) given by equation 
(22). The stationary solution \<j>"{x,t)\'^ (dashed line) is also depicted in the figure. 
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Figure 5. The main graph, as the previous one, shows the evolution of \^p"{x,t)\'^ 
(solid line) at a later time t = 0.3 fs. Note that the main part of the wave reaches the 
stationary solution |<^f, (a;, t)\^ (dashed line) at this short time. The small Sommerfeld- 
type precursor can be observed near the relativistic cutoff at a; = 90.0 nm. At the 
inset we show that the precursor (solid line), is well described by equation (22) 
(dotted line). 



Figure 6. Time evolution of \t(>^(x,t)\^ for different values of the position: xi = 0.4 
nm, X2 = 0.6 nm and xg = 0.8 nm. Notice that the transient behavior leading to 
the stationary regime exhibits an oscillating pattern similar to the diffraction in time 
phenomenon. 



Figure 7. This graph illustrates \ip!^{x,t)\'^ as a function of time for a fixed value of 
the position xi = 3.0 nm, beyond 2xp. Notice that in this case, the diffraction-type 
pattern clearly disappears. 



Figure 8. The time evolution of \ip!^(x,t)\-^ in order to exhibit the main peak 
shift of the wave, for different values of the position: a; = 0.1 nm (solid line), 
x = 0.3 nm (dashed line) and x = 0.5 nm (dotted line). The corresponding peak 
positions are ps, p2 and pi, respectively. Despite the fact that the three curves 
fulfill relativistic causality, the wave front main peak exhibits an apparent violation 
of Einstein causality. 




2 




Distance (nm) 



Fig. 5 




Fig. 6 




Time (fs) 




Time (fs) 



